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Abstract 

Effective electrostatic interactions between colloidal particles, coated with polyelectrolyte brushes 
and suspended in an electrolyte solvent, are described via linear response theory. The inner cores 
of the macroions are modeled as hard spheres, the outer brushes as spherical shells of continuously 
distributed charge, the microions (counterions and salt ions) as point charges, and the solvent as a 
dielectric continuum. The multi-component mixture of macroions and microions is formally mapped 
onto an equivalent one-component suspension by integrating out from the partition function the 
microion degrees of freedom. Applying second-order perturbation theory and a random phase 
approximation, analytical expressions are derived for the effective pair interaction and a one-body 
volume energy, which is a natural by-product of the one-component reduction. The combination of 
an inner core and an outer shell, respectively impenetrable and penetrable to microions, allows the 
interactions between macroions to be tuned by varying the core diameter and brush thickness. In 
the limiting cases of vanishing core diameter and vanishing shell thickness, the interactions reduce 
to those derived previously for star polyelectrolytes and charged colloids, respectively. 
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I. INTRODUCTION 

Polyelectrolytes are ionizable polymers that dissolve in a polar solvent, such as water, 
through dissociation of counterions. Solutions of polyelectrolytes are complex mixtures of 
macroions and microions (counterions and salt ions) in which direct electrostatic interactions 
between macroions are screened by surrounding microions. Polyelectrolyte chains, grafted 
or adsorbed by one end to a surface at high concentration, form a dense brush that can 
significantly modify interactions between surfaces in solution. When attached to colloidal 
particles, e.g., latex particles in paints or casein micelles in milk [3j], polyelectrolyte brushes 



can stabilize colloidal suspensions by inhibiting flocculation j^,|5|. Biological polyelectrolytes 
(biopolymers), such as proteins in cell membranes, can modify intercellular and cell-surface 
interactions. 

Conformations and density profiles of polyelectrolyte (PE) brushes have been studied by 
a variety of experimental, theoretical, and simulation methods, including dynamic light scat- 
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small- an gle neutron scattering 



,3, 



transmission electron microscopy P , neu- 
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tron reflectomet ry Id . surface adsorption atomic force microscopy [12j , self-consistent 
field theory |l3, |l4 JJ, Q, H fl 
theory j^j, Monte Carlo simulation 



211, Poisson-Boltzmann 

nn 



scaling theory 

, and molecular dynamics simulation [23|,|2^]. Com- 
paratively few studies have focused on electrostatic interactions between PE brush-coated 
surfaces. Interactions between neutral surfaces - both planar and curved (spherical) - with 



grafted PE brushes have been modeled using scaling theory |2(j], while interactions between 







charged surfaces coated with oppositely- charged PEs have been investigated for planar 
and spherical (colloidal) surfaces j3] via Monte Carlo simulation and a variety of theoretical 
methods. 

While microscopic models that include chain and microion degrees of freedom provide the 
most realistic description of PE brushes, simulation of such explicit models for more than one 
or two brushes can be computationally demanding. The purpose of the present paper is to 
develop an alternative, coarse-grained theoretical approach, based on the concept of effective 
interactions, which may prove useful for predicting thermodynamic and other bulk properties 
of suspensions of PE brush-coated colloids. Modeling each brush as a spherical shell of 
continuously distributed charge, we adapt linear response theory, previously developed for 
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29|, to derive effective electrostatic interactions. The 
theory is based on mapping the mult i- component mixture onto an equivalent one-component 
system of "pseudo-macroions" by integrating out from the partition function the degrees of 
freedom of the microions. Within the theory, microions play three physically important 
roles: reducing (renormalizing) the bare charge on a macroion; screening direct Coulomb 
interactions between macroions; and generating a one-body volume energy. The volume 
energy - a natural by-product of the one-component reduction - contributes to the total free 
energy and can significantly influence thermodynamic behavior of deionized suspensions. 

Outlining the remainder of the paper, Sec. |H] defines the model suspension of PE brush- 
coated colloids; Sec. 11111 reviews the linear response theory; Sees. II V I and IVI present analytical 
and numerical results for counterion density profiles, effective pair interactions, and volume 
energies in bulk suspensions; and finally, Sec. IVII summarizes and concludes. 



charged colloids 
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and PEs 



II. MODEL 



The system of interest is modeled as a suspension of N m spherical, core-shell macroions 
of charge — Ze (valence Z), core radius a, and PE brush shell thickness I (outer radius 
R = a + /), and N c point counterions of charge ze in an electrolyte solvent in volume V 
at temperature T (see Fig. HJ. The core is assumed to be neutral, the macroion charge 
coming entirely from the PE shell. Assuming a symmetric electrolyte and equal salt and 
counterion valences, the electrolyte contains N s point salt ions of charge ze and N s of charge 
— ze. The microions thus number iV + = N c + N s positive and N_ = N s negative, for a 
total of Nfj, = N c + 2N S . Global charge neutrality in a bulk suspension constrains macroion 
and counterion numbers via ZN m = zN c . Number densities of macroions, counterions, and 
salt ions are denoted by n m , n c , and n s , respectively. Within the primitive model of ionic 
liquids jsoj ]. the solvent is treated as a dielectric continuum of dielectric constant e, which 
acts only to reduce the strength of Coulomb interactions between ions. 

In PE solutions, the counterions can be classified into four regions: (1) those within narrow 
tubes enclosing the PE chains, of radius comparable to the Bjerrum length, A# = e 2 / (ek B T); 
(2) those outside of the tubes but still closely associated with the chains; (3) those not 
closely associated with the chains, but still inside of the PE shells; and (4) those entirely 
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outside of the macroions. Counterions in regions (l)-(3) can be regarded as trapped by the 
macroions, while those in region (4) are free to move throughout the suspension. Within 
region (1), the counterions may be either condensed and immobilized on a chain or more 
loosely bound and free to move along a chain. These chain-localized (condensed or mobile) 
counterions tend to distribute uniformly along, and partially neutralize, the chains. In our 
model, counterions in regions (1) and (2) act to renormalize the bare macroion valence. 
The parameter Z thus should be physically interpreted as an effective macroion valence, 
generally much lower than the bare valence (number of ionizable monomers). From the 
Manning counterion condensation criterion [J , according to which the linear charge density 
of a PE chain saturates at ~ e/A#, we can expect the bare charge in an aqueous solution to 
be renormalized down by at least an order of magnitude. 

The local number density profiles of charged monomers in the PE brushes, p m0 n{r), and 
of counterions, p c {f), are modeled here as continuous, spherically symmetric distributions. 
Charge discreteness can be reasonably neglected if we ignore structure on length scales 
shorter than the minimum separation between charges. Spherical symmetry of charge dis- 
tributions can be assumed if intra-macroion chain-chain interactions, which favor isotropic 
distribution of chains, dominate over inter-macroion interactions, which favor anisotropy. 

The density profile of charged monomers depends on the conformations of chains in the 
PE shells. Electrostatic repulsion between charged monomers tends to radially stretch and 
stiffen PE chains. Indeed, neutron scatteri ng e xperiments || on diblock (neutral-charged) 
copolymer micelles, as well as simulations |24|, provide strong evidence that the arms of 
spherical PE brushes can exhibit rodlike behavior. Here we assume the ideal case of fully 
stretched chains of equal length - a porcupine conformation j^] - and model the charged 
monomer number density profile by 

0, r>R 

a<r<R (1) 



<iiilr 2 
0, r < a, 

where r is the radial distance from the macroion's center. The model thus neglects configu- 
rational entropy of the PE chains, although it does include the entropy of the microions. 
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III. THEORY 



For the model suspension denned above, our goal is to predict distributions of microions 
inside and outside of the PE brushes and effective interactions between macroions. Adapting 
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2a and 



the general response theory approach previously applied to charged colloids 
PE solutions [29|, we reduce the multi-component mixture to an equivalent one-component 
system governed by effective interactions, and approximate the effective one-component 
Hamiltonian via perturbation theory. To simplify notation, we initially ignore salt ions. 
The Hamiltonian then decomposes, quite generally, into three terms: 

H = H m {{R}) + H c ({r}) + H mc ({R}, {r}), (2) 

where {R} and {r} denote collective coordinates of macroions and counterions, respectively. 
The first term in Eq. (|2*|) . 

i N m 

H m = H hc + - ^2 v mm (\Ri — Rj|), (3) 

is the macroion Hamiltonian, which includes a hard-core contribution ifhc (kinetic energy 
and hard-core interactions), and an electrostatic contribution due to the bare Coulomb pair 
interaction potential 

Vmm{r) = (4) 

er 

at center-center separation r. The second term in Eq. (J2J), 

H C = K C + - £ Vccdn-Tjl), (5) 

is the Hamiltonian of the counterions with kinetic energy K c interacting via the Coulomb 
pair potential v cc (r) = z 2 e 2 /er. The third term in Eq. (J2J), 

H mc = ^2^2v mc (\R i -r j \), (6) 

i=l j=l 

is the macroion-counterion interaction, which also may be expressed in the form 



J dRp m (R) J drp c (r)t; mc (|R - r|), (7) 
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where p m (R) = ^(R- — Rj) and p c (r) = X^i<5( r — r j) are the macroion and coun- 
terion number density operators, respectively, and t> mc (|R — r|) is the macroion-counterion 
interaction potential (to be specified in Sec. II V|) . 

The mixture of macroions and counterions can be formally reduced to an equivalent 
one-component system by integrating out the counterion coordinates. Denoting traces over 
counterion and macroion coordinates by ( ) and ( ) , respectively, the canonical partition 
function can be expressed as 

Z = ((exp(-pH)) c ) m = (exp(-/3tf eff )> m) (8) 

where (3 = 1/ksT, H c g = H m + F c is the effective one-component Hamiltonian, and 

F c = - k B T\n(exp\-(3(H c + H mc )}) (9) 



is the free energy of a nonuniform gas of counterions in the presence of the macroions. 

Now regarding tl 
perturbation theory 



Now regarding the macroions as an "external" potential for the counterions, we invoke 

Tin 
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23, |2J, |30j and write 

F c = F + f'dX (H mc ) x , (10) 



o 



where F = — /c s Tln (exp(— (3H c )) c is the reference free energy of the unperturbed coun- 
terions, the A-integral charges the macroions (i.e., the PE brushes) from neutral to fully 
charged, H mc represents the perturbing potential of the macroions acting on the counteri- 
ons, and (H mc ) x is the mean value of this potential in a suspension of macroions charged to 
a fraction A of their full charge. 

Two formal manipulations prove convenient. First, we convert the free energy of the 
unperturbed counterions to that of a classical one-component plasma (OCP) by adding and 
subtracting, on the right side of Eq. ([10jh the energy of a uniform compensating negative 
background j^jj, Eb = —N c n c v cc (0)/2. Here n c = N C /[V(1 — r^c)] is the average density 
of counterions in the free volume - i.e., the total volume reduced by the volume fraction 
Vtic = (47r/3)n m a 3 of the macroion hard cores - and t> cc (0) is the k — > limit of the Fourier 
transform of v cc (r). Equation (fTUj) then becomes 



F c = F OCP + I 1 dA (H me ) x - E b , (11) 

J 
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where Focp = Fq + Ej, is the free energy of a homogeneous OCP excluded from the colloidal 
hard cores. Second, we express H mc in terms of Fourier components: 

1 1 

(H mc ) x = — v mc (k)p m {k) (p c (-k)) A + — Jim [v mc (k)p m {k) (p c (-k)) A ] , (12) 
V k ^ V fc-o 

where v mc (k) is the Fourier transform of the macroion-counterion interaction and where 
p m (k) = Y,f=i exp(— ik • Rj) and p c (k) = J2f=i exp(— ik • Tj) are Fourier components of the 
macroion and counterion number density operators. The k = term is singled out because 
the number of counterions, N c = p c (0), does not respond to the macroion charge, but rather 
is fixed by the constraint of global charge neutrality. 

Further progress requires approximations for the counterion free energy. Applying second- 
order perturbation (linear response) theory, the counterions are assumed to respond linearly 
to the macroion external potential: 

p c {v) = J dr' X (r - r') J dr" p m (r> mc (r' - r") (13) 

or 

p c (k) = x{k)v m c{k)p m {k), MO, (14) 

where x{k) is the linear response function of the OCP. 

Combining Eqs. flT| - fl"4|) . the effective Hamiltonian can be expressed in the form of the 
Hamiltonian of a one-component pairwise-interacting system: 

H eS = H hc + - J2 i>eff(|Ri-R;|) + ^o, (15) 

where v e ff(r) = v mm (r)+Vi n d(r) is an effective electrosatic macroion- macroion pair interaction 
that augments the bare macroion interaction v mm (r) by a counterion-induced interaction 

Vmd(k) = x{k) [vmc(k)] 2 . (16) 

The final term in Eq. (fTK|) is the volume energy, 

7 1 

, (17) 



E = F OC p + ^limf ind (r) + N m lim 

2 r— >0 fc^O 



1 Z 

--n m v in d(k) + n c v mc (k) + —n c v cc (k) 



which emerges naturally from the one-component reduction. Although independent of the 
macroion coordinates, the volume energy depends on the average macroion density and 



thus has the potential to significantly influence thermodynamics. Evidently, the effective 
interactions depend on the macroion structure through the specific form of the macroion- 
counterion interaction v mc in Eqs. (fTEj) and ()17j). 

The OCP linear response function, proportional to the corresponding static structure 
factor S(k), may be obtained from liquid-state theory j^]. In practice, the OCP is weakly 
correlated, with coupling parameter T = \b/cl c 1, where a c = (3/47rn c ) 1 ' 3 is the coun- 
terion sphere radius. For example, for hard-sphere macroions of radius a = 50 nm, valence 
Z = 500, and volume fraction 7]y ic = 0.01, in water at room t emp erature (A# = 0.714 nm), we 
find T ~ 0.02. As in previous work on charged colloids j^H^J and polyelectrolytes j^J, we 
adopt the random phase approximation (RPA), which is valid for weakly-coupled plasmas. 
The RPA equates the OCP two-particle direct correlation function to its exact asymptotic 
limit: c^(r) = — j3v cc (r). Using the Ornstein-Zernike relation, S(k) = 1/[1 — n c c^(k)}, the 
linear response function then takes the form 

X (k) = -(3n c S(k) = - (1+ ^ c /jfc2) , (18) 

where k = y/4nn c z 2 XB is the Debye screening constant (inverse screening length). Note that 
the screening constant, which involves the density of counterions in the free volume, naturally 
incorporates the excluded volume of the macroion cores. With x(k) specified, the counterion 
density can be calculated from the macroion-counterion interaction and Eq. (JTljl for a given 
macroion distribution (see Sec. lIVI) . Finally, salt ions can be easily incorporated by introduc- 
ing additional response functions |28j . The pair interaction and volume energy are then mod- 



ified only through a redefinition of the Debye screening constant: k = y47r(n c + 2n s )z 2 \B, 
where n s = N s /[V(l — rjh c )] is the average number density of salt ion pairs in the free volume. 

Generalization of response theory to incorporate leading-order nonlinear microion re- 
sponse entails three-body effective interactions, as well as corrections to the effective pair 
potential and volume energy (3^]. Nonlinear effects are generally significant, however, only in 
concentrated, deionized suspensions of highly charged macroions 32j and are here ignored. 
It has been shown that response theory, combined with the RPA, is formally equivalent to 



Poisson-Boltzmann theory 32]. Both approaches rely on mean- field approximations that ne- 
glect microion fluctuations and predict microion distributions of the same general form, aside 
from a distinction in the screening constant, which response theory corrects for excluded vol- 
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ume of the macroion cores. Advantages of response theory over Poisson-Boltzmann theory 



are its predictions of (1) the entire effective Hamiltonian, including the one-body vo 
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28, 



33 



34 



ume en- 
35,33, 



ergy, which is essential for a complete description of phase behavior 
and (2) a more accurate expression for the Debye screening constant that incorporates the 
macroion excluded-volume correction. 



IV. ANALYTICAL RESULTS 



For our porcupine model of a spherical PE brush with 1/r 2 monomer density profile 
Gauss's law gives the electric field around a macroion as 

Zer — a 



E(r) 



0, 



r > R 
a < r < R 
r < a. 



(19) 



Integration over r yields the electrostatic potential energy between a brush and a counterion: 

Zze 2 

(20) 



v mc (r) 



er > 



Zze' 



Zze z 



1 - 



In 



a 



Hi) 



R 



r > R 
a < r < R 
r < a, 

aecause v mc (r) is not uniquely defined inside 
3q , we choose a below by requiring that the 



where a is an arbitrary constant, which arises 
the hard core. Following van Roij and Hansen 
counterion density vanish inside the hard core. Fourier transforming Eq. (J2*U|) yields 

AnZze 2 



skH 



-G(ka, kR; a), 



where the function G(ka, kR; a) is defined as 



G(xi, X2] a) = sinc(x2) — sinc(xi) — a[x\ cos(xi) — sin(xi)], 



(21) 



(22) 



with sinc(x) = f$ du sin(u)/u. We can now calculate, in the dilute limit, the counterion 
number density profile around a single macroion, taking p m (k) = 1. From Eqs. ()14j) . 1)180. 
and 1)21)1 . the Fourier component of the density profile is 



Pc{k) 



K 



z 



z kl(k 2 + k 2 



-G(ka, kR; a), 



(23) 



Pc{r) 



Z K 

z Airlr 



which in real space takes the form 

S(/ca, kR; a) e~ Kr , r > R 

S(«a, nr; a) e~ Kr + Ec(— nr, — KR)smh(nr), a<r<R (24) 

[Ec(— Ka, —kR) + a(l + na)e~ Ka ] sinh(Kr), r < a. 

For simplicity, we have introduced two functions, S(xi,X2',a) and Ec(xi,X2), which are 
defined, respectively, as 



S(xi,X2] a) = shi(x2) — shi(xi) — a[x± cosh(xi) — sinh(xi)] 



(25) 



and 



Ec(xi, x 2 ) = Ei(x 2 ) - Ei(xi), (26) 

where shi(a;) = Jq du sinh(u)/w denotes the hyperbolic sine integral function and Ei(x) = 
/^oo due u /u is the exponential integral function. Now setting the counterion density to zero 
within the hard core [i.e., p c (r) = 0,r < a, in Eq. 1)240] fixes the constant a: 



a 



-Ec(— Ka, —kR). 



(27) 



1 + Ka 

Integrating Eq. ()24j) over the spherical shell volume of a PE brush yields the fraction of 
counterions inside a brush: 

1 + kR kR 



fn 



1 - 



kI 



S(ko, k.R; a) 



28 



From this expression, it is clear that the counterion distribution, within the model, is deter- 
mined entirely by two independent dimensionless parameters, na and kI, i.e., the ratios of 
the macroion core radius and brush thickness, respectively, to the Debye screening length. 
From Eqs. (JTfij) and f)21j). the induced electrostatic pair interaction is given by 



Vind(k) 

whose Fourier transform is 

find(r) = 



AnZ 2 e 2 



l 2 k\k 2 + K 2 



[G(ka,kR;a)f 



2Z 2 e 2 n 2 /■<*> sin(Ar) 2 



For nonoverlapping brushes, Eq. ()30j) can be reduced to the analytical form 



2^2 



Vind(f) 



Z 2 e 



er 



Z 2 e 



2„2 



S(na, kR; a) 
kI 



1 2 



r > 2R. 



(29) 



(30) 



(31) 
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After adding to Eq. (j31j) the bare Coulomb potential between the spherical macroions 
[Eq. (HJ)], the residual effective pair interaction is 



Z 2 e 



2^2 



v e fi(r 



S(na, kR; a) 
kI 



r > 2R. 



(32) 



Thus, within the coarse-grained PE brush model and at the level of linear response theory, 
nonoverlapping PE brushes are predicted to interact via an effective Yukawa pair potential 



of the same screened-Coulomb form as the long-range limit of the DLVO potential 



37J for 



charged colloids. This result is consistent with previous linear response results for charged 



hard spheres 



27 



which interact via the DLVO effective pair potential 



Z 2 e 



2^2 



kR 



and for PE stars 



which interact via 



r > 2R, 



(33) 



Z 2 e 



2^2 



shi(/ti?) 
kR 



r > 2R. 



(34) 



Note that the screening constant, k, in the pair potential depends on the total density of 
microions - inside and outside of the brushes - since all microions respond to the macroion 



charge. We do not consider here over 
between chains also should be included 



apping brushes, in which case steric interactions 



24 1 . 



Finally, the volume energy is obtained from Eqs. (fTTj) . (|2"Tj) . (|2T)|) . and (jSHJ), as 

Z 2 e 2 K 2 oo [G(ka,kR;a)] 2 „ T AT ,k B T 



En 



'OCP 



N„ 



dk ■ 

reel 2 Jo k 2 (k 2 + K 2 ) y ' ^ "~'2 

Assuming weakly-coupled microion plasmas, the OCP free energy is well approximated by 
its ideal-gas limit: 



(35) 



Focp = iV+K"+A+)-l] + iV„[ln(n_A 3 _)-l], 



(36) 



where A± are the thermal de Broglie wavelengths of the positive/negative microions. The 
physical interpretation of the volume energy is straightforward. The first term on the right 
side of Eq. (J35)) represents the entropy of free microions, the second term the electrostatic 
energy of microion-macroion interactions, and the third term accounts for the background 
substraction. If the macroion valence Z is allowed to vary with concentration (e.g., through 
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counterion condensation), then Eq should be supplemented by the macroion self energy. 
We emphasize that, because of its dependence on the average macroion concentration, the 
volume energy has the potential to influence thermodynamic phase behavior. As a check 
of the present results, it can be shown that in the two limiting cases of vanishing PE shell 

lard core diameter (a — > 0) 



thickness (I — > 0, with Z fixed) and, independentl y, v anishing 
all analytical results reduce to those given in refs. 23, |28j and 



V. NUMERICAL RESULTS AND DISCUSSION 



respectively. 



To illustrate applications of the theory developed above, we present numerical results for 
the case of monovalent counterions (z = 1) in aqueous suspensions at room temperature 
(Xb = 0.714 nm). Figure 121 shows the predicted counterion profiles around three different 
types of macroion, all of the same outer radius R = 50 nm, valence Z = 500, and reduced 
number density n m R? = 0.01, for a salt-free suspension. The chosen valence is within the 
upper limit suggested by charge renormalization theory [38] for this size of macroion: Z < 
O(10)R/Xb- For a star macroion, the counterion density diverges logarithmically towards the 
center 29|, while for brush-coated and bare hard-sphere macroions the counterion densities 
remain finite. Figure H3 displays the corresponding internal counterion fraction, i.e., fractional 
counterion penetration, as a function of kR. For fixed ratio of hard-core radius to outer 
radius, a/R, the internal counterion fraction increases monotonically with k,R, reflecting 
increasing permeability of the macroions to counterions with decreasing screening length 
(e.g., increasing salt concentration). On the other hand, when kR is fixed, the counterion 
penetration decreases upon thinning of the PE brush (increasing a/R). In the limit of 
vanishing brush thickness (l/R — > 0, a/R — ► 1), Eq. (J28J) reduces to 

Ul^0) = ^-Kl + O(f). (37) 
na + 1 

Counterions are predicted to penetrate PE brush-coated macroions less efficiently than stars. 

Penetration of macroions by counterions can strongly influence screening of bare Coulomb 
interactions. Thus, effective pair interactions between brush-coated macroions depend sensi- 
tively on the thickness of the PE brush. To illustrate, Fig.|3]shows the effective pair potential 
for the same three macroion types as in Figs. El and El and for two salt concentrations, c s = 
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M and c s = 100 fiM, corresponding to different Debye screening constants k. For identical 
system parameters, the strength of the Yukawa pair interaction for nonoverlapping brush 
macroions is intermediate between that for hard-sphere and star macroions. Figure com- 
pares the dependence of the macroion-size-dependent amplitude of v e g(r), r > 2R, on Debye 
screening constant for the three macroion types. The amplitude increases with kR for fixed 
ratio of hard-core to outer radius, while for fixed kR the amplitude increases from the star 
limit to the hard-sphere limit as the PE brush thins to infinitesimal thickness (a/R — *• 1). 

VI. CONCLUSIONS 

Summarizing, polyelectrolyte-coated colloids provide a valuable conceptual bridge be- 
tween charged colloids and polyelectrolytes. In this paper, linear response theory is applied 
to bulk suspensions of spherical colloidal particles coated with PE brushes. Assuming stiff, 
radially stretched PE chains, we model each brush as a spherically symmetric shell of con- 
tinuously distributed charge, the charge density varying with radial distance r as 1/r 2 . By 
formally integrating out the microion degrees of freedom, the Hamiltonian of the macroion- 
microion mixture is mapped onto the effective Hamiltonian of an equivalent one-component 
system. Predictions of the theory include microion density profiles, effective electrostatic 
interactions between pairs of (nonoverlapping) macroions, and a state-dependent one-body 
volume energy, which contributes to the total free energy. The theory presented here may 
provide a practical guide for choosing system parameters to achieve desired interactions. 

The main conclusions of this study are: (1) Trapping of counterions inside a spherical PE 
brush is highly sensitive to variations in the core radius, brush thickness, and Debye screening 
length of the solution. For fixed ratio of core to outer radius, the fraction of trapped counteri- 
ons increases monotonically with increasing outer radius or decreasing screening length. For 
fixed ratio of outer radius to screening length, the fraction of trapped counterions decreases 
monotonically from a maximum in the limit of vanishing core radius (PE star macroion) to 
zero in the limit of vanishing shell thickness (hard-sphere macroion). (2) Within the linear 
response approximation, the effective pair interaction between nonoverlapping macroions 
has a Yukawa (screened-Coulomb) form. (3) By varying core radius and brush thickness, 
effective interactions between PE brush-coated macroions can be tuned - in both amplitude 
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and range - between interactions for hard-sphere and star macroions. For fixed ratio of core 
radius to outer radius, the amplitude of the pair interaction increases monotonically with 
increasing outer radius or decreasing screening length, while for fixed ratio of outer radius 
to screening length, the amplitude increases monotonically from the star-limit to the hard- 
sphere limit. The range of the pair interaction, governed by the Debye screening length, 
depends on the hard-core volume fraction and so can be varied by adjusting the core radius. 

The range of validity of the coarse-grained model and linearized theory studied here, and 
the accuracy of the predicted Yukawa form of effective pair interaction, including amplitude 
and range, could be directly tested by future simulations of more explicit models of PE- 
grafted colloids. Our purely electrostatic model can be augmented by chain elasticity and 
entropy - essential for describing overlapping PE shells j^jj. The mean- field linear response 



theory can be refined to incorporate nonlinear microion response 



and microion correla- 



tions, beyond the random phase approximation. The theory also can be easily adapted to 
other macroion types, such as core-shell microgels 2^], Future work will explore ther- 
modynamic phase behavior, which we anticipate to be quite rich and tunable between that 
of charge-stabilized colloidal suspensions and polyelectrolyte solutions. 
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FIG. 1: (a) Polyelectrolyte (PE) brush-coated colloidal sphere and, (b) model considered here, in 
which the PE monomer charge distribution is assumed continuous and varying as 1/r 2 , a < r < a+l. 




FIG. 2: Counterion number density profiles of three types of spherical macroion of outer radius 
R = 50 nm, valence Z = 500, and reduced number density n m B? = 0.01 in water at room 
temperature (\b = 0.714 nm): PE brush-coated macroion [solid curve from Eq. (|24|)] . PE star 
[dashed curve from Eq. (20) of ref. and charged hard sphere [dot-dashed curve from Eq. (32) 
of ref. 2?|]]. For the brush-coated macroion, the hard-core radius is a = 25 nm and the PE shell 
thickness is I = 25 nm. 
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FIG. 3: Fraction of counterions [from Eq. I)28jl ] trapped inside PE brush as a function of the 
dimensionless parameter kR (ratio of outer radius to Debye screening length) for several values of 
a/R (ratio of core radius to outer radius). 
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FIG. 4: Effective electrostatic interactions between pairs of nonoverlapping macroions of outer 
radius R = 50 nm, valence Z = 500, and reduced number density n m R? = 0.01 in room temperature 
water (Xb = 0.714 nm) at salt concentrations (a) c s = mol/1 (kR ~ 0.95) and (b) c s = 100 /umol/1 
{kR ~ 1.9): PE brush-coated spherical macroions [solid curves from Eq. PE stars [dashed 

curves from Eq. (|34|)]: and charged hard spheres [dot-dashed curves from Eq. (|33[)]. 
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FIG. 5: Amplitude of Yukawa effective electrostatic interactions between pairs of nonover lapping 
brushes, stars and charged hard spheres vs. Debye screening constant, normalized to unity at 
kR = [from Eqs. (|52|) - (fM|) ]. 
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